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Abstract 

A nanomechanical resonator coupled to a tunnel junction is studied. The oscillator modulates 
the transmission of the junction, changing the current and the noise spectrum. The influence of the 
oscillator on the noise spectrum of the junction is investigated, and the noise spectrum is obtained 
for arbitrary frequencies, temperatures and bias voltages. We find that the noise spectrum consists 
of a noise floor and a peaked structure with peaks at zero frequency, the oscillator frequency and 
twice the oscillator frequency. The influence of the oscillator vanishes if the bias voltage of the 
junction is lower than the oscillator frequency. We demonstrate that the peak at the oscillator 
frequency can be used to determine the oscillator occupation number, showing that the current 
noise in the junction functions as a thermometer for the oscillator. 

PACS numbers: 73.23.-b, 85.85. +j, 72.70.+m, 03.65.Ta 
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I. INTRODUCTION 



The recent years have seen the shrinking of mechanical components to micrometer and 

further to nanometer size, spawning the new field of nanomechanics.-^ Manufacturing tech- 

ft 

niques borrowed from semiconductor chip manufacturing (see e.g. Ref. y|) as well as bottom 
up approaches utilizing nanotubes^ make it possible to produce nanomechanical resonators 
with resonance frequencies presently demonstrated up to 1 GHz. Nanomechanical resonators 
are of interest for several reasons. In proposed tests of the limits of quantum mechanics 
one would like to investigate the decoherence behavior of superpositions of macroscopically 
distinct states, of e.g. a nanomechanical resonator.-^ Typical nanomechanical resonators 
contain a macroscopic number of atoms, ~ 10 7 , making the amplitude of the basic mode a 
macroscopic observable. The signature of such superpositions is strongest for low oscillator 
occupation numbers. Presently obtainable resonator frequencies are high enough to make 
cooling to the ground state feasible. Experimental efforts are under way to reach the first 
goal necessary for these measurements, cooling a nanomechanical resonator to the ground 
state.^^ii 

Other applications include the use of nanomechanical resonators as ultra-sensitive force 
detectors. Mass detection with zeptogram resolution utilizing nanomechanical resonators 
has been realized only recently— Nanosized cantilevers can also be used to detect magnetic 
forces. Detection of a single electron spin using a nanomechanical cantilever has already 
been demonstrated.— Nanomechanical resonators can also find an application in the context 
of quantum computing. Coherent mechanical oscillators are suggested as coupling elements 
between phase qubits in a solid state quantum computer.— All the mentioned applications 
not only require the fabrication of a suitable oscillator but also a way to detect the motion 

of a nanomechanical resonator. Different schemes for detection have been proposed (see 

I 

e.g. Ref. Mh. The most promising candidates for sensitive readout are electrical devices, 
such as tunnel junctions or single electron transistors, incorporated on the same chip as the 
nanomechanical resonator.- 

In the light of the possible applications it is necessary to obtain a theoretical understand- 
ing of nanomechanical resonators interacting with electrical devices on a chip. Theoretical 
descriptions of charge dynamics influenced by an oscillator have, up to now, mainly used 
a master equation technique. Mozyrsky and Martin investigated the model, where the 
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transmission coefficient of a tunnel junction depends on the position of a nearby harmonic 
oscillator, in the zero temperature limit. They found that the oscillator acquires an effective 
temperature proportional to the junction bias voltage and also find the influence of the os- 
cillator on the junction current at zero temperature.— Clerk and Girvin calculated the noise 
induced by a harmonic oscillator in a tunnel junction for dc and ac bias at zero tempera- 
ture using a Markovian master equation.— The authors, together with Khomitsky found the 
current in a tunnel junction influenced by an oscillator for arbitrary system parameters.— 
The current and the noise power spectrum for an asymmetric junction in the high voltage 
limit were also calculated. Smirnov, Mourokh and Horing^ analyzed the position fluctua- 
tions in the stationary state of a nanomechanical oscillator coupled to a tunnel junction for 
an exponential dependence of the tunneling amplitude on the oscillator position as well as 
the current through the tunnel junction. The theoretical descriptions of a nanomechanical 
oscillator interacting with a single electron transistor concentrated on the effects the single 
electron transistor, acting as a non-equilibrium environment, has on the oscillator. Ro- 
drigues and Armour— derived a master equation for an oscillator-single electron transistor 
system and investigated stationary state properties and dynamics. Blencowe, Imbers and 
Armour^ as well as Clerk and Bennett^ considered the interaction of a superconducting 
single electron transistor with a nanomechanical resonator and discovered that for a particu- 
lar source-drain voltage the single electron transistor can cool the oscillator and investigated 
a regime where the damping constant becomes negative . 

In this article we consider a tunnel junction coupled to a harmonic oscillator as a model 
for a nanomechanical resonator interacting with a detector and concentrate on the features 
of the noise power spectrum that the oscillator induces in the tunnel junction. We are 
interested in this model system for several reasons: The induced noise provides a means to 
determine the temperature of the oscillator. A shot noise thermometer, utilizing a tunnel 
junction, was demonstrated by L. Spietz et at— and was shown to work over a temperature 
range from 50 mK to 25 K. The noise induced by a nanomechanical oscillator in a SET 
has been successfully used to detect the temperature of a nanomechanical oscillator in two 
recent experiments,-^ in the region where the occupation number of the oscillator was large. 

We also want to investigate if there exists a signature of the oscillator in the noise power 
spectrum even if the oscillator is in its ground state and the voltage is insufficient to excite 
the oscillator. A recent treatment of a similar system, considering a spin instead of an 
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oscillator, claims a non-vanishing contribution to the noise power under these conditions.— 
A similar prediction could be made by blithely extending the result obtained by a Markovian 
master equation (e.g. from Ref . 1 1 7h into the region where the bias voltage is smaller than the 
oscillator frequency. Revisiting the junction/oscillator system utilizing a different technique 
will give us opportunity to investigate the question of the seemingly non- vanishing noise. 

Also there has been a recent theoretical discussion about which current-current correlator 
is detected in a noise experiment. Lesovik and Loosen,— Aguado and Kouwenhoven,— as 
well as Gavish et al— argue, that a passive detector, e.g. a LC oscillator at zero temperature 
or a two-level system, can only detect the positive frequency part of the Fourier transform of 
the unsymmetrized current-current correlator. The oscillator coupled to a tunnel junction 
gives us the opportunity to revisit this question in the context of a more complicated system 
than a mere tunnel junction. 

In this paper we will apply a Green's function technique to calculate the noise power 
spectrum of a tunnel junction coupled to an oscillator in the approximation of weak coupling, 
but for otherwise arbitrary parameters. The article is structured as follows: In section [Til 
we introduce the model Hamiltonian, and in section [TTT] we consider the stationary state of 
the oscillator using a Green's function technique. In section [TV] we calculate the average 
current through the junction as well as the unsymmetrized noise power spectrum. We 
consider application of the results in section [V], discussing noise thermometry. We present 
the conclusions in section |VT1 Details of the calculations are presented in appendices. 



II. OSCILLATOR INTERACTING WITH A TUNNEL JUNCTION 

Let us consider the situation of a nanomechanical resonator, modelled as harmonic os- 
cillator, interacting with a measuring device, modelled as a tunnel junction. The oscillator 
modulates the transmission amplitude of the junction thus changing the current and noise 
characteristics of the junction. The biased junction in turn acts as a non-equilibrium environ- 
ment for the oscillator, driving the oscillator from its initial state into a stationary thermal 
equilibrium state, albeit with a temperature different from the environment temperature of 
the tunnel junction. The Hamiltonian of the model system is 

H = Ho + Ht + Hr + Hr (2.1) 



4 



where H is the Hamiltonian for the isolated harmonic oscillator with bare frequency Qb and 
mass m. A hat marks operators acting on the oscillator degree of freedom. The Hamiltonians 
H[ r specify the isolated left and right electrodes of the junction 

H^^eicld , H r = J2 £ rclc r (2.2) 

1 r 

where l,r label the quantum numbers of the single particle energy eigenstates in the left and 
right electrodes, respectively, with corresponding energies £i jr and annihilation and creation 
operators. The operator describes the tunnelling, 

H T = f + f^ , T = J2TiAc r (2.3) 

with the tunneling amplitudes, T lr = T r \, depending on the oscillator degree of freedom. 
Due to the interaction of the tunnel junction and the oscillator, the tunnelling amplitudes 
and thereby the conductance of the tunnel junction depend on the state of the oscillator. 
In the following we assume linear coupling between the oscillator position and the tunnel 
junction 

T\ r = V\r + W\ r X (2.4) 

where v\ r = v* x is the unperturbed tunneling amplitude and w\ r = w*j its derivative with 
respect to the position of the oscillator. 

To discuss the current and noise in the tunnel junction, the current operator is needed 

I = i(f-f^. (2.5) 

The tunneling Hamiltonian consists of a part independent of the state of the oscillator 
and a part that depends on the state of the oscillator, so that the tunneling Hamiltonian, 
Eq. (12.31) . can be presented on the form 

Ht = h v + xh w . (2.6) 

For notational convenience, we have introduced the symbolic notation 

h u = T u + Tj, T u = y^ttfrtjcr, u = v,w, (2.7) 

lr 

where the symbol U\ r can take the values V\ r or w\ r . Similarly we can write the current 
operator, Eq. (|2.5I) . as 

I = 3v + xjv>, (2.8) 
5 



with 



Ju 



(%-%!), u = v,w. 



(2.9) 



When calculating current and noise in the tunnel junction the following combinations of 
the model parameters t>i r and u>i r appear 
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where / is the Fermi function. Here and in the following, the transmission matrix elements 
are assumed real. 



III. STATIONARY STATE PROPERTIES OF THE OSCILLATOR 

In this section we consider the properties of the stationary state a harmonic oscillator 
reaches due to interaction with a tunnel junction. In the Keldysh technique (for a review 



see e.g. Ref. 



271 ). we introduce the contour ordered oscillator matrix Green's function 

D(t, r') = -i (T c (x h (t)x h (t'))) . (3.1) 



The subscript H refers to an operator in the Heisenberg picture. Each of the times r and 
t' belong to one of the two branches of the Keldysh contour from — oo to +oo, T c denotes 
the contour ordering operator that orders operators along the Keldysh contour. We will use 
the symbol r to denote times on the contour, whereas t denotes real times. The branch 
index makes D a 2 x 2 matrix. The Keldysh matrix defined by Eq. (13. ip can be linearly 
transformed to the "triangular" form: 



D 




(3.2) 



where D R , D A , and D K are the retarded, advanced and Keldysh Green's functions, respec- 
tively. 

For a stationary state, the elements of the Keldysh matrix are functions of only the 
difference of real times t — t', and the Fourier transformed oscillator Green's function satisfies 
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the matrix Dyson equation 

(Dz\u>)-IL(u>))D(u) = i (3.3) 

where Dq (u) = [m{uj 2 — £1%)], being the bare oscillator frequency, and the self-energy 
(polarization operator) is a matrix of the form 

(u R U K \ 

n= . (3.4) 

Assuming weak interaction of the oscillator with the tunnel junction, the self-energy can 
be taken to lowest order. Calculations, details of which can be found in appendix [Aj give 
the following expression for the polarization operator, 



U{u) = -iG ww ' +R+ W {u)l, (3.5) 




where the conductance G ww is defined in Eq. (12.101) and the second term, the real part of 
the self-energy, is given by Eq. (|Alllk for to of the order of the oscillator frequency, R^ w (u) 
can be replaced by a constant R+ w (lu) ~ -R+^O).— The function 

Sv(uj) = — — — coth — ^- 1 coth — — — , (3.6) 

where T is the temperature of the junction and V = elf, U being the applied dc-voltage, is 
proportional to the well-known value of the power spectrum of current noise of the isolated 
junction, see e.g. Ref. l28l. Solving the Dyson equation, Eq. (13.31) . the retarded and advanced 
oscillator Green's functions become 

D R (fo)=nr l - , \ 2 — , D\u) = (D R {u))\ (3.7) 

where 7 e = — < 3U R (u)/2muj, the damping coefficient due to the coupling to the junction is 

G ww , . 

7e = W (3 - 8) 

and the renormalized oscillator frequency is 

fi 2 = fi|- 7e 2 + ^L(0). (3.9) 
For the Keldysh component we obtain 



D K {u) = (D R (co) - D A (uo)) . (3.10) 



Additionally to the environment provided by the coupling to the tunnel junction a 
nanomechanical oscillator is also subject to an intrinsic environment, e.g. phonons, act- 
ing as a heat bath and leading to damping. This additional heat bath, which we take 
to have the same temperature, T, as the junction, can be added phenomenologically, or 
expli citly b y adding the interaction with a bath of harmonic oscillators (as introduced in 



Refs. 



mil 



171). The total damping coefficient for the harmonic oscillator will then be the sum 



of the damping coefficients stemming from the tunnel junction, j e , and the heat bath, 70, 
giving a total damping coefficient 7 = 7 e + 7o- As a consequence of the additional heat bath 
the relation between the oscillator Green's functions, Eq. ( 13.101) , is modified according to 

D*M = - B») + ^coth^) , (3.11) 

\ 7 uj 7 IL J 

with the damping coefficient 7 e given in Eq. (13.71) . The relation for the Keldysh Green's 
function, Eq. (13.111) . is characteristic of the oscillator interacting with the environment which 
is in a non-equilibrium but steady state, and only in the absence of a bias voltage it reduces 
to the fluctuation-dissipation relation for an oscillator in thermal equilibrium at temperature 
T. However, since the coupling of the oscillator to the junction is weak, 7 -C max(V, Q, T), 
the oscillator spectral function is peaked at its frequency fl, and Eq. (|3.1ip can be written 
in the standard form 

D k (uj) = (D r (uj) - D a (uj)) coth (3.12) 
where the temperature T* characterizing the stationary state of the oscillator is given by 

,,1!"=^,^"). (3.13) 

2T* V7 Q 7 2T J 
The temperature of the oscillator, T* , depends not only on the environment temperature 
T, but also on the bias voltage of the junction V, as well as the relative coupling strengths 
7 e /7 and 70/7. For zero bias voltage the effective temperature reduces to the environment 
temperature. A frequency dependent effective temperature in the context of nanomechanical 
systems was also discussed by Clerk,— as well as by Clerk and Bennett,— for coupling to a 
general environment. Clerk derived a Langevin equation for a harmonic oscillator coupled 
to a fluctuating force, and obtained a relation similar to Eq. ( 13.131) . defining a frequency 
dependent effective temperature. 
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IV. NOISE PROPERTIES OF THE JUNCTION 



In this section we shall study how a harmonic oscillator interacting with a tunnel junction 
influences the current noise in the tunnel junction in the steady state. The noise properties of 
the junction current can be obtained in perturbation theory. The noise spectrum is specified 
by the current-current correlation function 

(5i H (t)5i H (t')) = (i H (t)i H (t')) - 1 2 , (4.i) 

where the current operator is given by Eq. (12.81) and the subscript H refers to an operator 
in the Heisenberg picture. For calculational convenience we introduce the current-current 
correlator with the time arguments lying on the Keldysh contour 

s(T,T') = (r c (5i H (T)5i H (T f ))). (4.2) 

It can be written in the interaction picture as 

S(r,r') = (t c ( e - l ^^M/(r)/(r'))) - J 2 , (4.3) 

where c denotes the Keldysh contour. To obtain the noise spectrum from the current-current 
correlator Eq. (14.31) we introduce Keldysh indices i, j = 1,2, that label the contour, e.g., 
i = 1 for the forward contour or i = 2 for the backward contour, and revert to using the real 
times t and i! , so that 

S(t,t') -> S ij (t-t'). (4.4) 
Finally, taking the Fourier transform of Eq. (14.41) we obtain 

/oo 

dte™ (4.5) 
-oo 

For calculational purposes it is sufficient to consider the real-time unsymmetrized current- 
current correlator 

S<(t - t') = S 12 {t - t') = (I(t)I(t')) , (4.6) 

since all other correlators can be derived from it, e.g. S > (t — t') = (J(t')J(t)) = S < (t' — t). 

A. I-V characteristic 



First we calculate the average current / to second order in the tunneling amplitude 

I{t) = jdr (T c (H T {T)i(t))) . (4.7) 
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Inserting the tunneling Hamiltonian Eq. (12.61) and the current operator Eq. (I2.8D we get two 
contributions to the current 

where one part is given by the standard result for the tunnel junction 

Ivv (*vv V- (4.9) 

with the conductance given by Eq. (12.1 0f> . 

The contribution induced by the coupling to the oscillator, I ww , can be written as 

1 f°° 

Iww = ^-.G ww l duJ(u). (4.10) 



2i 

where the conductance is given by Eq. (|2 . lOj) and 



J(u) = VD K (u) - [D R (u) - D a (lu)] A v (u), (4.11) 



with the oscillator Green's functions found in section IHIl and 



A/x V + u ,V + uj V-uj V-uj 

AvM = coth — — coth (4.12) 

Since the oscillator Green's functions are peaked at the oscillator frequency the integral in 
Eq. (14.1 Op can be evaluated and using the expression for the Keldysh Green's functions, 
Eq. (13121) . we get 

I ww = \G WW [V + N* + V. (N* + 1)] , (4.13) 

where 

h 

Cr W w G wv] (4.14) 

and we introduced the short notation 

V± = V ±A v (n) (4.15) 
and the occupation number of the oscillator is given by the Bose function 

N * = ^T~ V ( 4 - 16 ) 

where the temperature of the oscillator, T*, is specified in Eq. ( 13.131 ). From Eq. ( 14.131) it can 
be seen that two oscillator-assisted tunnelling processes contribute to the current: tunneling 
electrons gaining energy from the oscillator and loosing energy to the oscillator, respectively. 
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B. Current-current correlator 



Next we turn to calculate the current-current correlator. For an isolated junction, the 
noise is given by the second order expression in the tunneling amplitude since the fourth 
order correction only introduces a small featureless correction. However, when the junction 
is coupled to a quantum system exhibiting resonant behavior, as in the case of the oscillator, 
the interaction of the system with the tunnel junction can markedly increase the fourth order 
noise correction at the resonance and combinational frequencies. We will be interested in the 
resonant contributions and will therefore consider the current-current correlator to fourth 
order. 

Expanding the expression for the current-current correlator, Eq. 04.31) , to fourth order in 
the tunneling amplitude we obtain 

S(t,t')=S {2 )(t,t') + S w (t,t'), (4.17) 

with the correlator to second order in the tunneling amplitude 

^(2)(r,r') = (T c (/(r)/(r'))) (4.18) 

and the correlator to fourth order in the tunneling amplitude 

S(4) (t,t') = 

- \ Jdn Jdr 2 (t c (^(n)^^)/^)/^')) ) l , (4.19) 

where the subscript L denotes that only linked diagrams contribute to the expression, since 
the disconnected diagrams are cancelled by the current squared term. We first investigate 
the second order in tunneling contribution to the current-current correlator, resulting in the 
noise floor. 



1. Noise floor 

In this section we are going to discuss the noise floor. The main contribution to the floor 
comes from the second order contribution to the noise spectrum. 

To second order in the tunneling amplitudes the current-current correlator 

S (2) (r, t') = S vv (r, t') + S ww (r, r'), (4.20) 
11 



consists of the correlator for the isolated junction 

S vv (r,r') = (T c {j v (r)j v (r'))) = -iU+ v (r,r'), (4.21) 
and a contribution induced by the coupling to the oscillator 

S ww (t, t') = (T c {j w {r)j w {r'))) (T c (x(r)x (r'))> , (4.22) 

and therefore 

S ww (r, r') = -U+ w (t, t')D(t, t'), (4.23) 

where and 11+^ are given by Eq. (1 Alh . The correlator for the isolated junction can be 
presented as 

S<{u)=G vv [S v (w)-u]. (4.24) 

This result has been previously obtained by Aguado and Kouwenhoven for the quantum 
point contact.— 

The symmetrized noise spectrum of an isolated junction 

S«(u) = S<(u) + S>(u), (4-25) 

becomes the well known result^ 

S£(w) =2G vv S v (uj), (4.26) 
where we used the property, S < (u) = S > (—u). In the limit of zero voltage the noise floor 



reduces to Johnson 
Ref. 



16 and Ref. 
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NTyquist noise. The Markovian master equation approach employed in 
is not able to reproduce the correct frequency dependence, but only 
captures the zero frequency noise, and gives a constant noise floor with the magnitude of 
S K (0). 

At low temperatures T < V, the noise spectrum S< v is controlled by the voltage V: at 
positive frequencies to > \V\ the noise is exponentially small, S^ v (uj) ~ 0, and S< v increases 
linearly with the distance from the threshold uj = \V\: 

3<{u) = G vv [(V - lu)6(V -u) + (-V - lu)6(-V - u)] . (4.27) 

Eq. (14.271) shows that at T = 0, the noise power S< v (uj) is proportional to the phase space 
volume available for tunnelling events, i. e., the total number of electron-hole states, with an 
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electron on one side of the junction and a hole on the other side, with the excitation energy 
oj. Obviously, this result holds only in the frequency range where the electron density of 
states can be considered as a constant. 

The contribution to the noise induced by the oscillator, Eq. (14. 231) . becomes according 
Eq. (j335 

/oo 
duo 1 [S v {uo l )-uj l }D < {uo-uo l ) (4.28) 
-oo 

where D < (uj) = [D k (lj)—D r (lj) + D(uj)}/2. Recalling that the oscillator Green's functions 
are peaked at the oscillator frequency, the resulting contribution to the noise can be presented 
as 

^(w) = ^G ww {N* [Sv(cu~) — cu-] 

+ (N* + l) [S v {u + ) -u + ]} (4.29) 

where u>± = u ± Q and the conductance G ww is given by Eq. 04.141) . 

The second order contribution of the oscillator to the noise is similar to that of an isolated 
junction: with a different coupling, G ww in the place of G vv , the noise is given by the 
same expressions apart from the frequency shift ±Q. The two terms in Eq. (j4.29p give 
the noise contribution due to processes of electron tunneling accompanied by absorption or 
emission of oscillator quanta, respectively. When the oscillator approaches the ground state, 
N* — > 0, and at low junction temperatures, T <C Q, V, the noise S< w (u) is seen to vanish 
for frequencies larger than V — Vt. 

As expected, there will be no contribution to the noise Eq. f)4.29p for positive frequencies 
at zero temperature and voltage |V| < since the voltage is insufficient to excite the 
oscillator. 

To summarize, the second order contribution to the noise spectrum consists of the well 
known noise spectrum of an isolated junction, and a part that depends on the state of the 
oscillator. The second order contribution to the noise spectrum is frequency dependent and 
changes on the scale of max (V, Q, T). 

2. Resonant contribution to the current- current correlator 

In this section we consider the resonant contribution to the noise spectrum, that is, sharp 
peaks in the noise power with a width given by the oscillator damping 7. Technically, the 
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j i j i 

Figure 1: The four diagrams proportional to v 2 w 2 that contribute to the peak at the oscillator fre- 
quency. The wiggly lines represent the oscillator propagator and the bubbles are defined according 
to Fig. [3 The short arrows labelled by ui indicate the external frequency, i, j, k and I are Keldysh 
indices. 

resonant contribution originates from the fourth order terms Eq. (14.191) . Among various 
Feynman diagrams generated by the expression in Eq. (I4.19p . we are therefore interested 
only in those that produce resonant features. The diagrams of interest are shown in Figs. 

where the wiggly lines represent oscillator Green's functions, while the bubbles labelled 
with (-) are antisymmetric combinations of bubbles of electron Green's functions as specified 
in Fig. El 

The diagrams with a single oscillator line in Fig. Q] generate sharp features in the noise 
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u — U\ 



Figure 2: The two diagrams proportional w 4 that contribute to the peaks at u = and uj = 20. 
The wiggly lines represent the oscillator propagator and the bubbles are defined according to Fig. 03 
The short arrows labelled by u represent the external frequency, the indices i, j, k and I are Keldysh 
indices. Integration over the internal frequency, uj\ is implicit. 




Figure 3: The diagrammatic representation of the junction Green's functions, Eq. (lAip . Each 
bubble represents a linear combination of two electron bubble diagrams, with the Keldysh indices 
i and j. Integration over the internal frequency, U\, is implicit. 

spectrum at uj = ±fi, while the two-line diagrams in Fig. [2] are responsible for noise peaks 
in the vicinity of u = and uj = ±2Q. 

An example of a 4-th order diagram whose contribution is a featureless function of fre- 
quency uj is shown in Fig.[H In this diagram, the frequency of the oscillator line is integrated 
over and therefore the resonant contribution is absent. Being small compared with the second 
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Figure 4: A diagram representing a non-bubble contribution to the noise spectrum. The solid 
lines represent left and right side electron Green's function, the wiggly line represents the oscillator 
Green's functions. Integration over oj\ and u>2 is implied. 

order contribution to noise, this diagrams and other diagrams of this type can be discarded. 

We will denote the resonant contribution of the fourth order diagrams to the current 
correlator by S< es - It is given by the sum 

SZ.(u) = S$ v r{«>) + S&(u), (4-30) 

where Sf£ 2 (u) and S^ 4 (u>) are the contributions of the diagrams in Fig. [Q and Fig. El 
respectively. 

Evaluating the expression represented in Fig. Q] by performing the summation over the 
contour indices k and I, and inserting the bubbles and Green's functions we get (see Appendix 
[B]for details) 

S< 2w *H = ~2iGl w VJ(u;), (4.31) 

where the function J{oS) is specified in Eq. (14.111) and the conductance given in Eq. (12.101) . 
This contribution to the noise spectrum is peaked at u = ±Q due to the resonant behavior of 
the harmonic oscillator spectral function. The dependence of the peak height on the system 
parameters will be discussed in section IIV CI 

The contribution containing two oscillator Green's functions in Fig. [2] is calculated in 
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Appendix iBl to give 

s~t1 roc 

S< 4 (cu) = / duj.JWJiu-uJ 

^ J -oo 
r~i2 roc 

+ ^ / duiAMAfa - wi) 

2 J-oo 

x [y - Ay(^i)] [V - A y (cj - CJi)] (4.32) 

where we introduced the oscillator spectral function A(w) = i [D R (u) — D a {uj)\ . The de- 
pendence of the contribution to the noise spectrum, Eq. (14.321) . on system parameters will 
be discussed in section IIV CI 



C. Resonant contributions to the noise spectrum 

The resonant contribution to the noise consists of peaks at zero frequency, uj = 0, the 
oscillator frequency, uj = ±Q, and twice the oscillator frequency, uj = ±20. Their heights 
depend inversely on the damping coefficient of the oscillator, 7, and their width is propor- 
tional to the damping coefficient. In this section we will discuss the properties of the peaked 
structure and its dependence on environment parameters like bias voltage and temperature. 
The peaked contribution to the noise can be written as 

S<» = S< 2 (u) + S^u) + S<(u) + S<(u) + S<(u), (4.33) 

where Sq(oj) describes the peaked contribution to the noise spectrum at zero frequency, 
S±i(w) the contributions at the oscillator frequency and S± 2 (uj) the contributions at twice 
the oscillator frequency. 

The peaks at the frequencies uj = ±Q originate from the bubble contribution, 55 2(^)5 
given in Eq. 04.311) . Inserting the explicit form of the oscillator Green's functions, Eqs. (13.71 
I3.12p . in Eq. (14.311) . and using that the oscillator Green's functions are peaked at the oscilla- 
tor frequency, we obtain for the resonant contribution to the noise spectrum at the oscillator 
frequency 

S±iH = 7 ^2— ^ p n- (4-34) 

The peak height at the resonance frequency is 

Pn = -& UW V [N*V + + (N* + 1) VC] , (4.35) 

7 
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where V± is given by Eq. (I4.15P and N* is the occupation number of the oscillator. The 
peak height scales inversely with the damping coefficient 7. For large voltages V ^> O the 
peak height is linear in the occupation number N* and quadratic in the voltage. The result, 
Eq. (14.351) . extends the result from the Markovian master equation calculation in Refs.QQ 
to voltages smaller than the oscillator frequency. 

The peaks at frequencies lu = and 00 = ±2f2 come from the fourth order contribu- 
tion quadratic in the oscillator Green's functions, Eq. (|4.32p . The remaining integration in 
Eq. (I4.32p can be done and the dominating contribution for weak damping comes from the 
poles of the oscillator Green's functions. Collecting the contributions to the peak at zero 
frequency we get 

S » = -^^0 (4.36) 

where the peak height P is 

^0 = ^G 2 WW V [N* 2 V+ + (N* + 1) V_] . (4.37) 

The peak height at zero frequency also scales inversely with the damping coefficient 7 and 
depends quadratically on the conductance G ww . The result for the peak at zero frequency, 
Eq. (14.361) , coincides with the result obtained by using a Markovian master equation approach 



as in Ref. 



171 . which therefore correctly captures the low frequency noise. 



The contributions at double the oscillator frequency differ for positive and negative fre- 
quency. In the vicinity of uj = 2Q, 

s ^ ^ Fl^ fai (4 ' 38) 

with the peak height given by 

P2n = G 2 WW ^VN* [V + N* + V. (N* + 2)] (4.39) 
47 



and for u ~ — 2Q 



with the peak height 



S - M= i„ + $ + W P - m (440) 



-an = Gl w i-V(N* + 1 



1 

47 

x [V+ (N* - 1) + V- (N* + 1)] . (4.41) 
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The peaks at double the oscillator frequency also show resonant behavior, i.e. the peak height 
increases with decreasing damping. For large voltages and occupation number, N* ^> 1, it 
depends quadratically on the oscillator occupation number and the voltage. 

The peak heights depend on the bias V, the environment temperature T and the relative 
coupling strengths 7 e /7 and 70/7. In the following we discuss the peak heights in some 
limiting cases. 



1. Dominant coupling to the tunnel junction 

If the coupling between the oscillator and the tunnel junction is much stronger than the 
coupling to the thermal environment, 7 e 3> 70, the occupation number of the oscillator is, 
according to Eq. (13. 13|) . given by 

a. For high temperatures T>7,0 the occupation number of the oscillator becomes 
N* « T/Q. The functions V± = V in this limit. The peak heights of the resonant contribu- 
tions to the noise spectrum are then at the oscillator frequency 

Pn = G 2 VW -V 2 ^ (4.43) 
7 \l 

and is linear in the environment temperature. The peak height at zero frequency is 



Po = ^WjQ > ( 4 - 44 ) 

and depends quadratically on the environment temperature, whereas the peak heights at 
twice the oscillator frequency become 

P±2n = G 2 WW ^V 2 (£) , (4.45) 

which also scales quadratically in the environment temperature. For dominant coupling to 
the junction 7 oc G ww and therefore all the peak heights becomes linear functions of the 
conductance. The peak heights for high voltages V ^> T, Q can be obtained from the peak 
heights at high temperatures by replacing T — > V/2 in Eqs. (14.431 l4~44l I4.45H . 
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b. At low temperatures, T <C f2, and low voltages V < Q, the occupation number 
approaches zero, N* 0. The peaks at the oscillator frequency as well as the peak at zero 
frequency and the peak at u> = 2Q disappear 

P = Pn = P2n = 0. (4.46) 
At uj — — 2Q we obtain a dip in the noise spectrum with depth 



P-2Q — -G 2 WW — V 2 . 



1 



Since in nanomechanical systems the coupling strength between the junction and the oscil- 
lator can be tuned, we will also discuss the case where the thermal environment dominates 
over the junction. 



2. Dominant coupling to the thermal environment 

Another limiting case to consider is the situation when the coupling to the thermal 
environment dominates over the coupling to the tunnel junction, 70 ^> 7 e . In the high 
temperature limit, T V, Q, we obtain the same results for the peak heights as in the 
case of dominant coupling to the junction, since for both cases the junction and the thermal 
environment act as thermal equilibrium environments with temperature T. For high voltages 
and temperatures much larger than the oscillator frequency, V ^> T ^> Q, we get 

j\T = ^M + ^»l. (4.47) 

The peak height at the oscillator frequency is 

,2 4 T/ 2 ( 7e \V\ , 70^ 



The peak height at zero frequency is 



The peak height at twice the oscillator frequency 



1 t,2 ( 7e I V" I , 7 T 



P ±2n ^-^ W + --) ,4.50) 



2 



20 



is also quadratic in the occupation number. All peak heights depend inversely on the damp- 
ing coefficient 7, since in the parenthesis 7 appears only in the relative coupling strengths 
7 e /7 and 70/7. 

For low temperatures T <C V, ft we distinguish two regimes. For low voltages, V < ft, we 
obtain the same results as in the case for strong coupling to the junction since N* ps 0. For 
voltages V 3> ft 



N* 



le\V\_ 

7 2ft' 

The functions V± ~ F ± ft and we obtain for the peak height at the oscillator frequency 



(4.51) 



G -V u 

vw „ 

7 



2|2iM 
7 2ft 



(4.52) 



for the peak height at zero frequency 



G 2 —V 2 

ww 27 



7 2ft y 



7 2ft 



(4.53) 



and the peak heights at twice the oscillator frequency 

1 



P2Q. — G,, 



'2 7 



-V 1 



7 2ft 



Je\V\ 

7 2ft 



(4.54) 



and 



P-2Q - G ww — V 



+ 1 



V 



l_e\V\ 

7 2ft 



-ft 



(4.55) 



JeM 

7 2ft 

We have obtained the peak heights of the resonant peaks in the noise spectrum of a tunnel 
junction coupled to a harmonic oscillator for arbitrary parameters V, T and 7 e ,7o- We find 
that if the oscillator approaches the ground state, iV* — > 0, the peaks at positive frequencies 
as well as the peak at zero vanish. In the next section we are going to discuss an application 
of the properties of the peaks in the noise spectrum of the junction: noise thermometry. 



V. NOISE THERMOMETRY 



In experiments trying to cool a nanomechanical oscillator to the ground state a diagnostic 
tool is needed to check the state of the oscillator and to confirm, eventually, that the oscillator 
really is in the ground state. Coupling the oscillator to an electrical device, e.g. a tunnel 
junction, gives a means to experimentally determine the state of the oscillator, since the 
oscillator influences the current and noise in the junction. In a noise thermometry setup the 
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noise that the oscillator induces in the tunnel junction is used to determine the temperature 
or occupation number of the oscillator. 

As demonstrated, the noise spectrum of a tunnel junction consists of three peaks, one 
at zero frequency, one at the oscillator frequency and one at twice the oscillator frequency, 
and the peak heights depend on the oscillator occupation number. If the oscillator couples 
weakly to the junction the highest peak is at the oscillator frequency Q. The relation for the 



peak height, Eq. (I4.35p . can be used to determine the oscil 



ator occupation number in an 



experiment. The usual experimental procedure (see Refs. is to measure the peak height 
as a function of temperature. Fig. [5] shows the expected outcome of such a measurement, 
the dependence of the peak height at the oscillator frequency, obtained from Eq. ( I4.35p , as 
a function of the environment temperature T for different bias voltages V and dominant 
coupling to the junction r y e /l = 10. The solid line shows the peak height for V = 0.1Q. 
It decreases monotonically for decreasing temperature, and becomes exponentially small at 
low temperatures. The dashed and the dotted line show the peak height for voltages larger 
than the oscillator frequency (V = 2fl and V = 5Q respectively). 

The height of the peak P^ can be conveniently scaled to a dimensionless quantity N c , 

N c = 1 Pn , (5.1) 
AG 2 V 2 

where the combination j/G^ is a property of the device. For high occupation numbers 
(i.e. high temperatures or high voltage, max(V,T) Q), the dimensionless peak height N c 
gives the oscillator occupation number: N* « N c . The normalization constant, AG 2 W V 2 /j 
in Eq. (15.11 ), can be read off as the high temperature slope in a plot of Pq vs. T/VL. In the 
region of low occupation numbers we obtain the relation for the occupation number 

N* = N e -~. (5.2) 

where V_, Eq. (14.151) . is a known function of the bias and temperature. 

We note that for low temperatures and voltages the oscillator occupation number is not 
simply proportional to the peak height at the oscillator frequency. The relative occupation 
number N*/N c is shown in Fig. [6] for dominant coupling to the junction and different bias 
voltages. At high temperatures the relative occupation number approaches unity, for low 
temperatures and low voltages it departs from unity. 

Another possibility to get information on the oscillator occupation number is from the 
peak heights at uj = and lj = 2fl using Eqs. (14.371 14.39B . The combination of junction 
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5 T/Q 



Figure 5: The peak height of the noise at the oscillator frequency, Pn, as a function of the 
environment temperature T for different bias voltages and dominant coupling to the junction, 
7e/7o = 10, and G vw = 0.1. The full line shows the temperature dependence of the peak height 
for V = O.lfi, the dashed line for V = 2Vt and the dotted line for V = 50. As the environment 
temperature approaches zero, the peak height vanishes for voltages V < $7, but stays finite for 
voltages V > £1. The square dots mark different occupation numbers of the oscillator as indicated. 

parameters V 2 G 2 llw /^ in Eqs. (|4.37l I4"l39j) can be obtained, e.g, from the high temperature 
slope in a plot of P vs - T 2 /fl 2 . 

Since the peak at the oscillator frequency is not isolated, but sits on a noise floor, the 
question of observability of the peak arises. As a measure of observability we define a signal 
to noise ratio, the peak height relative to the floor 

and demand r > 0.1 for the peak to be observable. The peak height, Pn, as well as the 
magnitude of the noise floor, S VV (Q) + S WW (Q), depend on A^*, the junction bias, and the 
junction temperature as well as the conductances G vv , G vw and G ww . To investigate the 
observability of the peak, we use Eq. f!5.3p to find the occupation number A^* for a given 
bias voltage, junction temperature and a set of conductances, so that r = 0.2. The solid 
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1 2 3 4 5 T/Q 

Figure 6: The relative occupation number of the oscillator N* /N c as a function of the environment 
temperature T for different bias voltages V and dominant coupling to the junction, 7 e /7o = 10, 
with G vw = 0.1. The full line shows the relative occupation number for low voltage, V = 0.1O, the 
dashed line and the dotted line show the relative occupation number for high voltages, V = 2Q and 
V = 50, respectively. For high temperatures T 3> f2, the relative occupation number approaches 
unity. For low temperatures the relative occupation number can depart considerably from unity to 
approach 1/(1 + ^/Q/j e V). 

lines in Fig. show the occupation number necessary to fulfill r = 0.1 as a function of 
voltage and junction temperature. To detect an oscillator close to its ground state high bias 
voltages and / or low junction temperatures are necessary, as seen in Fig. On the other 
hand, high voltages heat the oscillator, so that a compromise between heating and readout 
has to be found. To illustrate the heating effect, the dash-dotted lines show the occupation 
number for dominant coupling to the junction, 'j e /lo — 100. If one would like to read out 
the oscillator occupation number when N* < 0.1, without introducing heating, voltage and 
temperature are confined to the shaded region on the left of Fig. This example shows 
that a good signal to noise ratio for the readout and a minimal heating of the oscillator are 
complimentary requirements. 
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Figure 7: The solid lines show are a contour plot of the minimum occupation number necessary to 
observe the peak at the oscillator frequency with a S/N ratio r > 0.2, as a function of voltage and 
temperature with relative conductances G vv /G vw = 10 and G ww /G vw = 0.1 . The dashed lines are 
a contour plot of the occupation number of the oscillator for dominant coupling to the junction, 
7e/7o = 100. The shaded region to the left marks the parameter regime that allows a readout of 
the oscillator occupation number without heating the oscillator. 

VI. CONCLUSION 

We have considered a nanomechanical resonator interacting with a dc-biased tunnel junc- 
tion. We model the resonator as a harmonic oscillator and the interaction is introduced via 
the modulation of the tunneling amplitude by coupling to the harmonic oscillator position. 
Employing a Green's function technique we calculated the properties of the stationary state 
of the oscillator, obtaining that the coupling to the junction introduces damping and heats 
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the oscillator. The expression Eq. (13.131) . for the temperature of a harmonic oscillator cou- 
pled to a tunnel junction was also obtained in Ref. [L7j using a Markovian master equation 
approach. The reason for the coinciding results derived by two different techniques is due 
to the same weak coupling approximation made in both derivations. Both, the perturbative 
Green's function calculation presented here and the Markovian master equation can only 
be applied for weak coupling between the environment and the oscillator. In the master 
equation approach the weak coupling leads to a separation of time scales of the environment 
evolution and the evolution of the density matrix, resulting in a Markovian master equation. 
Here the same argument is used in the frequency domain: The width of the oscillator spectral 
function is small on the scale of the frequency dependence of the environment correlation 
functions. 

The current- voltage characteristic of the junction has been obtained. The stationary state 
current is seen to consist of two contributions. One contribution is the current through an 
isolated junction, whereas the other contribution, dependent on the state of the oscillator, 
stems from oscillator assisted tunneling. The expression obtained for the dc-current is in 
accordance with the result previously derived using a master equation technique.— 1 ^ The 
additional current arising from the influence of the oscillator on the junction transmission, 
Eq. (14.131) . vanishes if the oscillator is in the ground state and the voltage across the junction 
is smaller than the oscillator frequency. We observe therefore that the zero point fluctuations 
of the oscillator position do not affect the electric current. 

The main part of the paper presents the calculation of the electric noise due to oscillator- 
assisted tunnelling. The unsymmetrized current-current correlator has been evaluated for 
arbitrary frequencies, bias voltages, and environment temperatures. The noise spectrum 
consists of a smooth noise floor and a peaked resonant structure. 

As a function of frequency u, the noise floor varies on the scale max(T, V, O). The 
contribution of the oscillator is shown to vanish at positive frequencies when the oscillator 
is in the ground state and the bias voltage is smaller than the oscillator frequency. The 
expressions Eqs. ( 14.241 14.29D for the second order noise generalize the result obtained by a 



Markovian master equation derived in Ref. 



171 . In the Markovian master equation calculation 



a constant noise floor was obtained, whereas the Green's function calculation presented here 
captures the frequency dependent noise floor. As expected, in the low frequency limit, we 
recover the result obtained by the Markovian master equation approach. 
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The resonant structure in the noise spectrum, S < (u), consists of peaks at zero frequency, 
uj = 0, the oscillator frequency, u = ±fi, and twice the oscillator frequency, uj = ±20. Being 
present only for a finite voltage across the junction, the resonances are of non-equilibrium 
origin, and their intensities at positive and negative frequencies are not related by the 
detailed balance relation.— The peaks at the oscillator frequency have the same height 
for positive and negative frequencies, whereas the peaks at double the oscillator frequency 
are asymmetric. The peaks at u — ±f2 stem from processes involving a single vibrational 
quantum, and their height is therefore linear in the oscillator occupation number N*, whereas 
the processes leading to the peaks at u = and u = ±2Q involve two quanta and the peak 
heights are quadratic in N* . If the oscillator is in the ground state, the peaks at positive 
frequencies vanish in the limit where the bias voltage is smaller than the oscillator frequency 
but a dip in the noise spectrum remains for negative frequencies. 

For bias voltages smaller than the oscillator frequency and the oscillator in the ground 
state, we find no oscillator dependent noise at positive frequencies: Both the oscillator de- 
pendent contribution to the noise floor as well as the peaks at positive frequencies disappear 
from the unsymmetrized noise spectrum, S < (u). The oscillator contribution to the current- 
current correlator at negative frequencies remains finite even for voltages smaller than the 
oscillator frequency. 

To understand this result we consider a setup for measuring current fluctuations. Lesovik 
and Loosen^ as well as Gavish, Levinson and Imry,— introduced a damped harmonic oscilla- 
tor with resonance frequency flo as a meter of current fluctuations. For linear coupling of the 
oscillator to the current in the junction, they obtained the deviation of the meter position 
fluctuations, (X 2 ), from the equilibrium value in terms of the current-current correlators of 
the junction, S < and S > (see Eq. (I4.6P ). and the meter occupation number, Nq , 

5 (X 2 ) = A [(JVno + 1) S<(n ) + N Qo S>{tto)} , (6.1) 

where A is the coupling strength and Nq q is the Bose function with temperature T. Note 
that the argument of the noise spectrum S < ' > , the oscillator frequency Q is positive. Let 
us consider the meter fluctuations in different temperature ranges. For large detector tem- 
peratures T>n , and consequently Nq q ^> 1, the meter fluctuations are proportional to 
the symmetrized current-current correlator of the junction, S > + S < . A passive detector 
- a detector at low temperature T f2 , when Nq C 1 - measures the unsymmetrized 
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current-current correlator S< at positive frequency. In accordance with these arguments,— >2 6 - 
only the part of S < at positive frequencies describes physical noise, understood as random 
flow of energy from the system to environment. 

We can now understand our results for the current-current correlator S < (u) with the 
help of Eq. (16. ip . A passive detector detects only the positive frequency part of the current- 
current correlator. Our calculations show vanishing S < (u) at positive frequencies 10 in the 
case when the oscillator is in the ground state, and thus the ground state does not contribute 
to the physical noise in the tunnel junction in compliance with general expectations. 

In experiments with nanomechanical resonators utilizing electrical devices as detectors, 
noise properties can function as a diagnostic tool in determining the state of the resonator. 
We have shown how peaks in the noise power spectrum can act as a measure for the oscillator 
occupation number and discussed the criteria for observing these peaks. 

In this paper we presented a calculation of the noise induced by an oscillator in a tun- 
nel junction. The obtained results are valid at arbitrary parameters and thus also in the 
important region where the oscillator approaches the ground state. 

Appendix A: JUNCTION GREEN'S FUNCTIONS 

In this section we are going to calculate the junction Green's functions that appeared in 
the calculation of the properties of the oscillator stationary state, the average current and 
the noise. The junction Green's functions encountered are 



where the notation was introduced in section [III and the contour times discussed in section 
IHIl In terms of the tunneling operators T and T' , Eq. (12.31) . taken in the interaction picture, 
we have 



n~ u ,(Vi,T 2 ) 



i (T c (7i u (Ti)/v(-r 2 ))) , 



(T c (h u (Ti)j u '(T 2 ))) , 



(Al) 




(A2) 
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Inserting the tunneling operator from Eq. (12.71) we get the symmetric and antisymmetric 
combinations of left and right electrode electron Green's functions 



n* u /(ri, r 2 ) = -i ^ Ui r u' Tl [G x {t 2 - r^G^n - t 2 ) 



lr 



±G r (r 2 - r^G.in - r 2 )] 



(A3) 



where we introduced the electron Green's function for the right and left electrodes, for 
example 

Gi(ti,t 2 ) = -i (T c (^(rOcJCra))) . (A4) 
The Green's functions of interest can be obtained by standard techniques. For example, 

Kw±(tiM) = -i \(ri{t 2 )r u {t x )) ± (r^Xlih)) 



is given in terms of its Fourier transform 



lr 

x {[1 + n(e r - ei)] 5(e r - e x + V - w) 
± n(e r - ei)5(e r - ei + V + lu)} , 



(A5) 



where n{u) is the Bose function. For voltages V <C Ep we obtain the approximate relation 



Eir«ir<i [/(ei)-/(e r )]5(e r - ei + y) 
= ^£>r^r> (-^)5(e r - ei ), 



(A6) 



as tunneling amplitudes depend only weakly on their arguments, so that 



n< n , ± M = iG uu ,[{V-u)[l+n{V-u))] 



±{V + uj)n(V + 



where we introduced the conductance 



G uu > = 27T ^ UlrKl \~^T) ^ ~ ^ 
lr ' 



(A7) 
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The retarded, advanced and Keldysh bubbles 



n^±(ti-* 2 ) 
nf„'±(£i - *2) 



-ieih-h) [([^(ti),^(t 2 ) 
±<W(ti),r^)]>] 
-iefa-h) [([7;(ti),rj(t 2 ) 
±<W(ti),r tt ,(t 2 )]>] 
-< [({^(to,^^) 

±<W(*0,^(* 2 )}>] 



(A8) 



can be calculated in a similar way to give 



II* (u;) II* (u; 

o rrfH 



(A9) 



with 



C>) = -iGuwV + R- uu ,{oo) 
n„V(w) = -2zG^A y (a;) 



(A10) 



and the advanced part given by 11^ (u>) = [ir^(u;)] , with the functions Sy(uj) and A v (uj) 
defined in Eqs. (13.61 14.121) respectively, and similarly U* u ,(uj) is given by Eq. (13 . 5f) . 
The real parts of the retarded H^ u ,(lu) are given by 



R± u ,(u) = 2tt ]T u Ir <i [/(ei) - /((*)] (P- ± P + 



(All) 



lr 



with P± = 1/ (ej — e r ± V + uj). To estimate the real parts, we ignore the weak energy 
dependence of the transmission amplitudes and the electron density of states are constants 
up to a cut-off frequency of the order of the Fermi energy Ep. For frequencies and voltages 
much smaller than the cutoff frequency, u, V <C Ep, we can approximate the reactive parts 
of the response functions Eq. (lAlip to give 



to 



Hip 



(A12) 



This concludes our calculation of the junction Green's functions. 
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Appendix B: FOURTH ORDER DIAGRAMS 



In this appendix we give the algebraic expressions for the fourth order diagrams con- 
tributing to the current-current correlator. Let us consider, for example, the contribution 
containing bubbles that is linear in the oscillator Green's function, i.e., the diagrams depicted 
in Fig. Q] 




kl 

+Tl n (0)Il k2 (0)D 12 (Lu) 
+U n (uo)U k2 (0)D k2 (uj) 

+Tl a (0)Tl k2 (uj)D kl (uj)) , (Bl) 

where the full oscillator Green's function appears, accounting for the influence of the inter- 
action with the tunnel junction and the matrix 



(B2) 



introduces the correct signs for forward and backward contour (to simplify notation, we 
suppress subscripts v and w on the right hand side). Performing the summation over the 
contour indices k and I, and inserting the bubbles and Green's functions we obtain Eq. (14.311) . 

Similarly, the contribution from diagrams containing two oscillator Green's functions, see 
Fig. El becomes 

■S^M = ^iA'V^i), (B3) 

kl J 

A w (o;,a;i) = U l l(ui)U k2 (oj 1 )D kl (u; 1 )D 21 (oj - u x ) 

+ U l l(u - wi)n* 2 (o;i)D 2, (a; - w^D*^)- 

Doing the summation and inserting the bubbles and oscillator Green's functions we obtain 
Eq. (I4.32D . where we neglected a non-resonant contribution that is small compared to the 
second order contributions to the noise floor, Eqs. (14.291 |4~24"1) . 

All diagrams not containing bubbles are of the non-resonant kind, i.e. oscillator Green's 
functions always appear under frequency integration. A typical non-bubble contribution to 
the noise spectrum is shown in Fig. [4] and given by 
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kl liri l 2 r 2 

(B4) 

Since there are no resonant peaks present, contributions of this type have to be compared to 
the second order contribution to the noise floor, Eqs. (14.241 14.29D and since they are always 
of higher order in tunneling they can be neglected. 
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171 is not present here, since 



The additional pumping contribution to the current found in Ref. 
for simplicity we take the junction to be symmetric implying real transmission matrix elements. 
32 We have not explicitly included terms in the self energy, that correspond to a shift in the 
equilibrium position of the oscillator. The shift can be absorbed by measuring the oscillator 
coordinate from the new equilibrium position. We have to keep in mind, though, that the 
coordinate independent part of the tunneling amplitude, v r \, has to be changed accordingly. 
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